ASTIA  (TISVV) 

ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


Department  of  Defense  contractors  must  be 
established  for  ASTIA  services,  or  have 
their  need^to^know  certified  by  the  cogni“ 
zant  military  agency  of  their  project  or 


contract. 


AEDC-TR-61.4 


HYPERSONIC  FLOW-BLAST  ANALOGY 


By 

J.  Lukasiewicz 
VKF,  ARO,  Inc. 


June  1961 

Program  Area  750A,  Project  8952,  Task  89512 
ARO  Projects  300116,  360105,  360106 

Contract  No.  AF  40(600)-800  S/A  11(60-110) 


AF  -  .AEDC 
Arnold  AFS  Tenn 


AEDC-TR-61-4 


ABSTRACT 


Two-dimensional  and  axisymmetric,  inviscid,  hyper¬ 
sonic  flows  about  simple  slender  bodies  are  considered  with 
particular  reference  to  pressure  distributions  and  shapes 
of  shock  waves.  Approximate  solutions,  based  on  blast 
analogy,  are  derived  and  compared  with  more  exact,  theo¬ 
retical  calculations  and  with  experimental  results.  The 
validity  of  correlation  parameters  predicted  by  the  blast 
analogy  and  the  range  of  their  applicability  are  investigated. 
Near-exact  theoretical  results  available  for  read  airflow  at 
Mach  numbers  from  15  to  19  are  compared  with  the  perfect 
gas  data. 
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NOMENCLATURE 


A,  B 
a 

a,  b 

C 

D 

d 


f,  g,h 

Jo 

M 

Mn 

Mr 

p 

P' 

o 

Psh 

% 

R 

Ro 

r 

S 

t 

u 

V 


X 


X 


Functions,  Eq.  (62) 

Velocity  of  sound 
Constants 
Function,  Eq.  (64) 

Drag  coefficient 
Drag 

Plate  thickness  or  cylinder  diameter 

Explosion  energy  per  unit  area  of  the  surface  of  shock 
front  when  R  =  1 

Functions,  Eq.  (3) 

Constant 

Mach  number 

Mach  number  component  normal  to  shock  front 
Shock  Mach  number  (  =  V/a^) 

Static  pressure;  pressure  on  plate  or  cylinder  surface 

Pitot  pressure 

Shoulder  pressure 

Pressure  at  shock 

Distance  from  shock  to  origin 

®a/P»c)  ^ 

Radial  or  normal  co-ordinate 
Distance  measured  along  shock 
Time 

Velocity  component  in  the  free-stream  flow  direction 
Velocity  of  blast  wave 

Velocity  component  normal  to  the  free-stream  flow 
direction  or  flow  velocity  caused  by  explosion 

Distance  from  nose  in  the  free-stream  flow  direction 

Distance  from  shoulder  in  the  free-stream  flow  direction 


9 


AEDC-TR-61.4 


y 

a 

y 

K 

P 

T 

x' 

SUBSCRIPTS 

oo 

SUPERSCRIPTS 

(  )' 

(  )" 


< 

Constant;  =  0  for  plane  flow,  =  1  for  axisymmetric  flow 
Functions,  Eq.  (13) 

Ratio  of  specific  heats 
r/R 

Density 

Body  thickness  parameter 
Functions,  Eq.  (8) 

Function,  Eq.  (76) 


Free  stream 


First  approximation 
Second  approximation 
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INTRODUCTION 


This  report  deals  with  the  problem  of  inviscid  hypersonic  flow  of 
air  over  simple  slender  bodies.  Two-dimensional  and  axisymmetric 
flows  are  considered  with  particular  reference  to  pressure  distributions 
and  shapes  of  shock  waves.  Approximate  solutions,  based  on  blast 
analogy,  are  derived  and  compared  with  more  exact,  theoretical  calcu¬ 
lations  and  with  experimental  results.  The  validity  of  correlation  param¬ 
eters  predicted  by  the  blast  analogy  and  the  range  of  their  applicability 
are  investigated.  Near-exact  theoretical  results  available  for  real  air¬ 
flow  at  Mach  numbers  from  15  to  19  are  compared  with  the  perfect  gas 
data. 

The  author  wishes  to  thank  Messrs.  H.  W.  Ridyard,  R.  E.  Geiger, 
and  the  General  Electric  Company  for  permission  to  use  their  unpublished 
flow  field  computations  and  Mr.  Vernon  Van  Hise  of  the  NASA  Langley 
Research  Center  for  making  his  original  data  available.  His  thanks  are 
also  due  to  Miss  P.  Mitchell  and  Mrs.  B.  Majors  for  their  help  with  the 
preparation  of  this  report. 


BASIS  OF  BLAST  ANALOGY 


The  equations  of  small  disturbance  theory  of  steady,  hypersonic  flow 
(Ref.  1)  may  be  written  as: 


Continuity: 

UooPx  +  +  apv/r  =  0 

Momentum: 


+  VVj,  +  Pj./p  =  0  V 


Energy: 


(1) 


with  a  =  0  for  plane  and  a  =  1  for  axisymmetric  flow. 


As  pointed  out  by  Van  Dyke  (Ref.  1)  and  Hayes  (Ref.  2),  these  equa¬ 
tions  do  not  involve  the  streamwise  velocity  u,  and  therefore  the  flow  in 
the  transverse  planes  can  be  treated  independently. 


Manuscript  released  by  author  April  1961. 
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The  equations  of  one-dimensional,  unsteady  flow  are: 
Continuity: 

Pf  +  (pv)j.  +  apv/r  =  0 

Momentum: 

''t  +  vv,  +  p^/p  =  0  L 

Energy: 

(pp“y)j  +  v(pp'^)j,  =  0 


(2) 


Hayes  (Ref.  2)  pointed  out  that  the  steady  flow,  as  described  by  the 
small  disturbance  theory,  is  analogous  to  the  unsteady  flow,  since,  by 
putting  =  dx/dt,  Eq.  (1)  becomes  identical  with  Eq.  (2).  In  other 
words,  for  an  observer  stationary  with  respect  to  the  undisturbed  fluid, 
the  small  disturbance  theory  flow  is  given  by  the  unsteady  solution  in 
one  fewer  space  co-ordinates,  so  that  flow  in  the  plane  normal  to  the 
direction  of  motion  can  be  treated  independently  of  x . 


If  a  solution  could  be  obtained  for  the  simpler  unsteady  problem, 
then  the  corresponding  steady-flow  problem  could  also  be  solved;  this 
procedure  has  been  generally  known  as  the  blast-wave  analogy  method. 
Before  this  method  is  applied  herein,  it  is  pertinent  to  recognize  the 
assumptions  inherent  in  the  steady  and  unsteady  solutions  involved  and 
to  examine  their  compatibility. 

The  hypersonic  small  disturbance  theory  is  applicable  to  slender 
bodies,  such  that  the  body  thickness  parameter  r<<  1  (since  terms  of 
the  order  of  are  neglected),  and  at  large  Mach  numbers,  such  that 
M^r  >  1  .  Since  the  longitudinal  velocity  component  is  eliminated  in  this 
approximation,  the  drag  of  a  body  is  given  entirely  by  the  energy  of  the 
transverse  flow:  the  drag  work  expended  by  moving  the  body  to  a  given 
transverse  plane  is  equal  to  the  increment  of  energy  of  the  transverse 
flow  up  to  that  same  plane.  The  energy  of  the  flow  in  the  transverse 
plane  per  unit  length  will  not  change  when  the  slope  of  the  body  surface 
is  zero;  therefore,  within  the  hypersonic,  small  disturbance  approxima¬ 
tion,  the  energy  of  the  transverse  flow  will  remain  constant  downstream 
of  the  forebody  or  nose  when  it  is  followed  by  a  zero-slope  afterbody, 
such  as  a  flat  plate  or  cylinder. 

Turning  now  to  the  case  of  unsteady  flow,  approximate  solutions  for 
the  spherical,  cylindrical,  and  plane  unsteady  flows  with  constant  energy 
have  been  obtained  for  cases  in  which  the  shock  wave  which  contains  the 
flow  field  is  strong,  the  shock  Mach  number  being  large. 
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Consider  a  configuration  consisting  of  a  blunt  nose  followed  by  a 
zero-slope,  flat  plate  or  cylindrical  afterbody  at  a  large  flight  Mach 
number.  For  such  a  body,  the  hypersonic  small  disturbance  theory 
will  not  apply  in  the  nose  region,  since  r  «  i  is  violated  in  that  region. 
Also,  it  will  not  apply  far  downstream,  where  the  shock  has  decayed 
and  flow  deflections  are  negligible,  because  no  longer  is  >  1  .  There¬ 
fore,  for  a  body  of  this  type,  the  hypersonic,  small  disturbance  theory 
will  apply  in  a  region  of  certain  limited  extent  downstream  of  the  nose. 

In  this  region,  the  body  surface  slope  is  zero,  the  energy  in  the  trans¬ 
verse  plane  is  constant,  and  hence  the  unsteady  solutions  are  applicable. 
The  success  of  the  blast-wave  analogy  will  thus  largely  depend  on  the 
accuracy  of  the  solutions  of  unsteady,  constant  energy  flows. 


APPROXIMATE  SOLUTIONS  OF  UNSTEADY,  PLANS  AND  CYLINDRICAL, 
CONSTANT  ENERGY  FLOW 


Following  the  work  of  Taylor  (Ref.  3),  Sakurai  (Ref.  4)  obtained  first 
and  second  approximations  to  unsteady,  constant  energy,  plane,  cylin¬ 
drical,  and  spherical  flows.  In  this  section,  Sakurai's  results  are  sum¬ 
marized  for  use  with  the  blast  analogy  method,  and,  in  general,  his 
notation  is  adhered  to. 


FIRST  APPROXIMATION 


Approximate  blast-wave  solutions  are  constructed  in  the  form  of 
power  series  in  y  =  where  Mjf  is  the  shock  Mach  number.  In  the 

first  approximation,  all  powers  of  y  except  zero  are  neglected;  that  is, 
a  strong  shock  is  assumed.  The  following  equations  result: 


v/V 


(o) 


f(/<) 


P/Poo 

p/Poo 


(o ) 

g(/<)/y 


y(R  /R) 


id  +  1 ) 


=  J„ 


(3) 


with  R^  =  (E^/p_^)  ^  where  E  is  the  explosion  energy  per  unit  area 

of  the  surface  of  the  shock  front  when  R  =  1 .  ,  and  h  are 

functions  of  k  =  r/R,  y  and  a  ,  whereas  is  a  constant  depending  on 
y  and  a  only. 
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They  are  given  by: 


(o) 

Kk)  =  Kly  +  Ax" 


o  / 

g(K)  = 


1  _  ^  (n  -  1) 


y  +  1 


2/^  +  (3a  +  l)y  -  (a  +  l) 
(1  -  a)y  +  3a  +  1 


^  ^  +  1.  / 

y  +  l  y-1  [y  +  l-  K 


(n  -  1) 


2(2a  +  1  +  y) 
(1  -  a)y  +  3a  +  1 


h(K) 


y  -  1 


with  A  =  i_  ^  ~  ^  .  n  =  (2  -  a)y"  +(3a  +I)y  -  1 


y  y  +  l 


y  -  1 


^  (4) 


The  variation  of  shock  position  and  pressure  with  time  are  obtained 
from  Eq.  (3)  with  V  =  dR/dt .  Thus 


ta, 

¥ 


^  =  ^/T„  (-5-)  " 

„  °  a  +  3  V  Ro  / 


3)/2 


(5) 


and 


R_ 


=  yfTo 


a  +  3 


(f)/fe) 


(a  +  3)  /  [2(a  +  1)] 


(6) 


The  boundary  or  shock  conditions  (k  =  i)  are  given,  to  this  approxi¬ 
mation,  by 

V  ^  2  P  1  ^  2y  P  ^  y+l 

V  y  +  1  ’  Poo  y  +  1  ’  y  -  1  (7) 

The  values  computed  by  Sakurai  (Ref.  4)  at  r  =  0  (/c  =  o)  for  y  =  1.4 
are  given  in  Table  1.  For  other  values  of  y,  is  given  in  Table  2,  and 
other  functions  can  be  computed  from  Eq.  (4). 
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SECOND  APPROXIMATION 

In  the  second  approximation,  all  powers  of  y  except  zero  and  one 
are  neglected,  and  the  following  equations  result: 

(o)  (o) 

v/V  =  Hk)  +  (k  -  Kk)  )<!'  (k)  y 
(o) 

P/Poo  =  (1  +  ^  (k)  y)/y 

p/Poo  =  h(K)  (1  +  yWy) 

y(R„/R)  (1  +  A^y) 

where  ,  if ,  and  x  are  functions  of  k  ,  y  ,  a  ,  and  \  depends  on  y  and  a  . 
The  shock  position  is  again  obtained  in  terms  of  time  as 

R/Hq 

=  ^/o  [(R/Ro)“^“  (9) 

and  the  pressure  ratio  in  terms  of  shock  location  as 

p/Poo  =  [  (Ro/R)  ^  -  -A)  ]  (10) 


For  R^oo  or  at  infinitely  large  time,  Eq.  (10)  gives  at  r  =  0  (k  =  0) 


P/Poo  =  -  "A) 

=  0.692  for  a  =  0 

=  0.606  for  a  ~  1 


y  -  1.4 


(11) 


For  plane  wave  (a  =  o),  Eq.  (9)  gives 

ta.»/Ro  =  (12) 

where 

/3  =-\Jif(  l  +  f)  -  logg  (Jl  +  i'+  j 

^  (13) 

and  V  =  -J^Ai  (R/R^  ) 


For  cylindrical  wave  (a  =  i). 


ta«,/Ro  =  I  1  -  [  1  -  JoA,  (R/R„)^  ]  j  (a.  yjX) 


(14) 
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In  this  approximation,  the  exact  conditions  at  the  shock  are  satis¬ 
fied  for  velocity  and  pressure,  that  is 


^  =  _ 2 _  (1  _  ) 

V  y  +  1 

P  -  1  y  -  1 


Poo  y  +  1  y  y  +  1 

the  above  equations  corresponding  to  Eq.  (8)  with 


and 


f(l)  =  _ ? _ 

y  +  1 


(o)  (o) 

3>(i)  =  -  f(i)/(i  -  f(i)) 


(15) 

(16) 


The  exact  density  ratio  across  the  shock  is 

2 


y  -  1 


(17) 


This  is  not  equivalent  to  the  corresponding  Eq.  (8),  except  for  small  y 
(large  ),  when  it  is  expanded  to 


(18) 


with 


and 


(o)  V  +  1 
h(l)  =  ^  ^ 


y  -  1 


yd)  = 


y  -  1 


The  coefficients  in  the  second  approximation  are  given  in  Table  1 
for  r  =  0  (k  =  0)  and  y  =  1.4. 


TABLE  1 

VALUES  OF  COEFFICIENTS  IN  THE  FIRST  AND  SECOND  APPROXIMATIONS  TO  THE 
BLAST  SOLUTION,  y  =  1.4,  k=0  (SAKURAI,  REFS.  4,  5) 


a 

(o) 

f(o) 

(o) 

g(o) 

(o) 

h(o) 

Jo 

O(o) 

^(o) 

y(o) 

A, 

0 

0 

0.  455 

0 

1.  696 

CO 

00 

CO 

-0.617 

3.  86 

-2.  138 

1 

0 

0.  424 

0 

0.  877 

-3.  5 

-0.56 

3.5 

-1.  989 
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TABLE  2 

VALUES  OF  Jo  (SAKURAI,  REF.  4) 


0 

1 

1.  2 

3.  024 

1.  547 

1.  3 

2.  147 

1.  102 

1.  667 

1.  137 

0.  585 

SOME  PROPERTIES  OF  APPROXIMATE  SOLUTIONS 

In  the  first  approximation,  the  ratio  of  pressures  at  distances  and 
is  constant  at  all  times  and  is  given  by  (see  Eq.  (3) ), 

Pk/P/<j  =  g('<i)  /g(k) 

For  y  =  1.4,  the  values  obtained  for  the  ratio  of  pressure  at  the  shock 
to  pressure  at  the  origin  are  given  in  Table  3. 


TABLE  3 

RATIO  OF  PRESSURE  AT  SHOCK 
TO  PRESSURE  AT  ORIGIN,  y  =  1.4 


a 

0 

1 

P  /p 

2.  56 

2.  68 

The  difference  in  the  pressure  ratio  p/p^,  as  given  by  the  first  (Eq.  (3) ) 
and  second  (Eq.  (8))  approximations,  is  constant  at  a  given  /<  and  y  and  is 
equal  to 

(o) 

'A('<)>0  (20) 

At  K  =  0  and  for  y  =  1.4,  this  amounts  to  0.  281  (a  =  0)  and  0.  237  {a  =  l). 


The  relationship  between  the  position  of  the  shock  and  the  shock  pres¬ 
sure  ratio  is  given  by 


(o)  P<x= 

g(l) 


(first  approximation) 


(21) 


17 


and 


(22) 


Jo 


(o) 

g(l) 


+  [a,  -  ^(1)]  gU)  I 

(second  approximation) 


For  y  =  1.4,  the  second  approximation  is 

a  =  0,  (R/R^)  (pjj/p^  -  2.33)  =  0.69 

a  =  1,  (R/R^)=  (Pll/Pc«  -  2.16)  =  1.33  J 
In  the  first  approximation,  the  numerical  term  in  the  bracket  is  omitted. 


(23) 


The  variation  of  pressure,  density,  and  velocity  with  distance  (k) 
from  the  origin  can  be  judged  from  Table  4,  in  which  values  of  coef¬ 
ficients  for  the  first  and  second  approximations  are  listed.  It  is  apparent 

that  pressure  {g^°\  ifr)  exhibits  negligible  variation  for  k<  o.4.  The  density 

variation  (h^°\  y)  is  also  small  near  the  origin,  whereas  velocity  4>) 
varies  nearly  linearly  with  k. 


TABLE  4 

VARIATION  OF  COEFFICIENTS  FOR  FIRST  AND  SECOND  APPROXIMATIONS  WITH  k, 
FOR  PLANE  AND  CYLINDRICAL  FLOW,  y  =  1.4 


a 

K 

0 

0.  2 

0.  4 

0.  6 

0.  8 

1.  0 

- - 

0 

0.  455 

0.455 

0.461 

— - - 

0.  496 

0.  631 

1.  167 

1 

0.  424 

0.  424 

0.425 

0.441 

0.  531 

1.  167 

0 

-0.  617 

-0. 573 

-0. 260 

0.  451 

-^0.  143 

1 

-0.  56 

- 

- 

-0.  45 

0.  092 

-0.  143 

J^(o) 

0 

0 

0.  039 

0.  225 

0.  669 

1.  766 

6.  000 

1 

0 

0.001 

0.  019 

0.  153 

0.  783 

6.  000 

X 

0 

3.  86 

3.  633 

3.  109 

1.  674 

-5. 000 

1 

3.  5 

“ 

- 

“ 

3.  002 

-5. 000 

^(o) 

0 

0 

0.  143 

0.  287 

0.  437 

0.  611 

0.  833 

1 

0 

0.  143 

0.  286 

0.  432 

0.597 

0.  833 

0 

-3.  86 

-3. 869 

-3.  803 

-3.  842 

-5. bob 

1 

-3.  5 

■" 

-3.542 

-5. 000 
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APPROXIMATE  STEADY-FLOW  SOLUTIONS  BASED  ON  BLAST  ANALOGY 


The  steady,  blast  analogy  approximations  are  obtained  from  the 
approximate  constant  energy,  blast  solutions  by  superposition  of  free- 
stream  velocity  =  x/t  and  by  expression  of  the  energy  of  unsteady 
flow  in  terms  of  body  wave  drag.  Thus  the  plane  blast  solution  is  anal¬ 
ogous  to  two-dimensional  flow  over  a  blunted  plate,  and  the  cylindrical 
blast  solution  corresponds  to  the  case  of  axisymmetric  flow  about  a 
blunted  cylinder.  In  each  instance,  first  and  second  blast  approxima¬ 
tions  are  used  to  obtain  approximate  steady-flow  solutions. 

In  obtaining  steady-flow  solutions,  the  shock  shape,  R/d,  and  the 
pressure  distribution,  Po/Po^-  on  the  body  surface  as  a  function  of  distance 
along  the  body,  x/d,  are  the  main  concern.  These  solutions  are  strictly 
applicable  only  to  flat  plates  and  cylinders  at  zero  incidence.  The  body 
pressure  p^,  denoted  below  by  p ,  is  taken  at  k  =0,  and  the  body  thick¬ 
ness  is  thus  neglected.  In  view  of  the  invariance  of  pressure  with  k  and 
the  small  density  near  the  origin  (see  Table  4)  this  procedure  seems 
reasonable. 

Static  pressure  and  other  parameters  of  flow  inthe  transverse  plane 
are  easily  obtained  from  blast  solutions,  provided  the  shock  shape  has 
been  determined.  Since 


Mr  =  (dR/dx)  (24) 

is  given  from  the  shock  slope,  distributions  of  transverse  pressure, 
density,  etc.  ,  follow  from  Eqs.  (3)  and  (8). 

It  should  be  noted  that,  although  the  second  blast-wave  approxima¬ 
tion  gives  essentially  exact  conditions  at  the  shock  surface  (see  Eqs.  (15) 
through  (18) ),  the  corresponding  blast  analogy  gives  correct  values  of 
pressure,  etc. ,  at  the  shock  only  when  the  component  of  Mach  number, 
normal  to  the  shock,  s. 


(dR/ds) 


(25) 


closely  approximates  This,  however,  is  consistent  with  the  assump 

tions  of  the  hypersonic  small  disturbance  theory. 


For  the  analogous  steady-flow  solutions,  the  blast  energy,  E^.  is 
expressed  in  terms  of  drag.  From  the  definition  of  E^, 


and 


E, 


=  -L  D„ 


"a=  1 


D 


27T 


a=  1 
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where 

=  npse  drag  of  plate  per  unit  length  of  leading  edge 

and 

=  nose  drag  of  cylinder 


In  terms  of  drag  coefficient,  C  , 


E 


,  ^  =  y  jL  p  Cn 

a  =  o  4  oo  oo  D^_q 


R 


‘’a  =  0 


y  4  Cd 


«  =  0 


and 


E 


i  ,  =  y  _dl  p  Cn 

a=l  ^  16  °°  ~  ‘^a=l 


R 


‘’a=l 


y  iL_  M 


(26) 


where  Cq  is  based  on  plate  thickness,  d,  or  cylinder  frontal  area,  jrdVd, 
and  free-stream  dynamic  pressure. 


FIRST  APPROXIMATION 


From  Eqs.  (5)  and  (6),  the  following  expressions  for  shock,  shock 
slope,  and  pressure  distribution  are  obtained. 

Plane  flow  (a  =  o); 


or 


d 


R/d 


0.774 


(x/ d)  ' 


(27) 

(28) 


(29) 


for  y  =  1.4; 


dR 

dx 


0.516 


r  1/3 


-1/3 


(30) 

(31) 
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and 


(o)  , 

g(o)  M!  C 


D 


=  0.121  Mt,  [CoAx/d)]'^’ 

for  y  =  1,4. 

From  Eqs.  (29)  and  (33), 

p/p^  =  0.0936  Cp/CR/d) 


for  y  =  1.4. 


Axisymmetric  flow  (a  =1): 

=  0.795  Cp*''"' (x/d) 
or 


for  y  =  1.4; 


and 


=  0.397  Cp  1/“  (x/d)-‘/2 


=  0.067 

oo 


y/(x/d) 


for  y  =  1.4. 

From  Eqs.  (37)  and  (41), 

p/p  =  0.0423  Cn  /(R/d)" 

for  y  =  1.4. 


(32) 

(33) 

(34) 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
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SECOND  APPROXIMATION 

From  Eqs.  (10)  to  (14),  the  following  expressions  for  shock  shape 
and  pressure  distribution  are  derived: 

Plane  flow  (a  =  0): 

In  this  case,  explicit  expressions  for  shock  shape  and  pressure 
ratio  were  not  obtained.  Thus,  for  shock  shape: 

t  -  Mi  (43) 

o 

=  0.066  P  Cn  M’  (44) 

for  y  =  1.4 

where  p  and  v  are  given  by  Eq.  (13)  and 

V  =  -(4J„A^/y)  (R/d)/  (M^  Cp)  (45) 

=  10.36  (R/d) /(M^  Cp) 

for  y  =  1.4. 

The  shock  slope  is  given  by 

(lit  1 

d  X  2 


=  0. 

From  Eq.  (10), 

^  =  -■^i/[(p/p„„)/g(o)  +  (K  -  ^)]  (49) 

=  0.973 /(p/p_^  -  0.692)  (50) 

for  y  =  1.4  which,  together  with  Eq.  (43)  gives  pressure  distribution. 

In  the  limit  when  (t,  R,  x)-»oo,  the  pressure  ratio  is  given  by  Eq,  (11) 
(=  0.692  for  y  =  1.4). 

The  function  p  (Eq.  (13)),  when  expressed  in  terms  of  pressure 
ratio  p/p^  (Eq,  (50) ),  can  usually  be  approximated  to  better  than  one  per¬ 
cent  by  an  expression  of  the  form 

/3  =  a(p/p^)^ 
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where  a  and  b  are  given,  for  y  =  1.4,  in  Table  5  for  various  ranges 
of  p/Poo. 


TABLE  5 

VALUES  OF  CONSTANTS  a  AND  b  FOR  VARIOUS 
RANGES  OF  p/p^ 


Range  of 

P/Poo 

a 

b 

2  -  4 

1.  30 

-1.  87 

4-20 

0.  87 

-1.  59 

20  -  100 

0.69 

-1.  52 

Using  this  form  of  fi ,  Equation  (43)  can  be  written 

(51) 

=  0.066a  Cp  C  (p/Poo)' 

(52) 

for  y  =  1.4; 

or 

Poo  L  ya  Cp  J 

(53) 

■  [  ] 

(54) 

for  y  =  1.4 

which  is  analogous  to  the  first  approximation,  Eq.  (32), 

A  more 
distribution 

convenient  form  of  the  second  approximation  for  pressure 
can  be  derived  as  follows: 

From  Eqs.  (32),  (43),  and  (49),  the  pressure  ratio  (p/p^)',  as  given 
by  the  first  approximation  at  a  given  x/d,  can  be  expressed  in  terms  of 
the  second  approximation  pressure  ratio  (p/p^)"  at  the  same  x/d  by 

[(P/Pj'  =  -(2/3)^/’  ]  ,/d 

(55) 

and 

(p/p.,c)'  =  0.743/1-^/’ 

(56) 

for  y  =  1.4 

and  is  therefore  a  function  of  (p/p„„)''  only. 
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The  values  of 

[(p/Poo)"  WPoc)']^/^, 

as  obtained  from  Eq.  (56),  are  tabulated  in  Table  6  and  plotted  in  Fig.  1 
in  terms  of  the  difference  (p/p^)"  -  (p/poc)'  vs  ip/p^)". 


TABLE  6 

PRESSURE  RATIOS  AT  SAME  x/d 
(a  =  0,  y  =  1.4,  FIRST  AND  SECOND  APPROXIMATIONS) 


(p/Poo)"x/d 

0.  692 

0.  8 

1 

2 

5 

10 

(p/Poo)  x/d 

0 

0.  191 

0.  428 

1.  473 

4.  493 

9.  503 

From  Fig.  1,  it  is  evident  that  the  simplified  expression 

(p/Pco)"  =  (p/Poo)'  +  0-56  (57) 

=  0.121  M:,  [Co/(x/d)]  \  0.56  (  58) 

gives  the  second  approximation  pressure  ratio  to  better  than  1.5  percent 
when  p/p^  2  1,  y  =  1-4. 


For  the  shock  shape  (Eqs.  (10)  and  (26)), 


n_ 

d 


p/Poo 

(o) 

_  g(o) 


+  (K 


(59) 


=  Co/ [  10.65 (p/p^)  -  7.37]  (60) 

for  y  =  1.4,  and  using  the  simplified  expression  (Eq.  (58) ), 

n/d  ^  _ 0.774 _ 

|^Co/(x/d)J  -1.09  (61) 

for  y  =  1.4. 

For  p/p^  >  1,  this  equation  gives  values  to  within  4.  5  percent  of  the  second 
approximation  (Eq.  (44)). 

Axisymmetric  flow  (a  =  1): 


J-  =  A^/X  /  1  -  2B  JL 
d  V  d  V  d 

with 

A  =  (yCo/jJ^A/V'^ 


(62) 
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and 

B  = 


hence 


for  y  =  1.4. 

^  =  A[C/2  -  B/C] 

dx 

with 

c  = 

^ _  (o)  - 

P  ^  'Ty  g(°)  VCp 

Poo  S’/jT' 


1  - 


2^1 /j7 
vT" 


x/ d 


Mio\^ 


(o) 

-  g(o)  (Ai 


-A) 


=  0.067 


X  /d 


x/ d 


+  0.606 


(63) 


(64) 


(65) 


(66) 


Using  the  same  procedure  as  in  the  case  of  plane  flow,  simpler 
equations  for  the  second  approximation  are  derived.  From  Eqs.  (40) 
and  (65), 


(o) 

_V=  ('_p_V 

1  _  ^1  g(o) 

/  \  Poo/ 

4(p/Pj„)/ 

“  (t)'  1 

'l  +  0.211/(p/p^ 

(o) 

g(o)  ~  Ip) 


.606 


(67) 

(68) 


for  y  =  1.4. 


Numerical  results  are  given  in  Table  7  and  Fig.  1.  From  Fig.  1,  it  is 
evident  that  the  simplified  expression 

(t)" 

=  0.067  +  0.44  (70) 


gives  the  second  approximation  pressure  ratio  to  better  than  one  percent 
when  p/p^  >  1,  y  =  1.4. 
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TABLE  7 

PRESSURE  RATIOS  AT  SAME  x/d 
(a  =  1,  y  =  1.4,  FIRST  AND  SECOND  APPROXIMATIONS) 


(p/Poc)" 

0.  606 

0.  958 

1.  432 

2.  415 

5.  404 

10.  4 

(p/p.)' 

0 

0.  5 

1 

2 

5 

10 

For  the  shock  shape,  Eqs.  (10)  and  (26)  give 

^  ^«>V^[p/poo  +  (^1  -  ’/')g(o)] 

=  0-206  M^y/c;;(p/p^  -  0.606)"^/^  (72) 

Using  Eq.  (69),  the  simplified  expression  is 

JL  =  0.795 

[m^  /H^/U/d)-  2.478]'/= 

For  p/p^  >  1,  this  expression  gives  values  to  better  than  one  percent  of 
the  second  approximation,  Eq.  (63). 


SUMMARY  OF  BLAST  ANALOGY  SOLUTIONS  FOR  y  =  1.4 


For  convenience,  the  above-derived  results  for  y  =  1.4  are  listed 
as  follows: 

Plane  Flow 


First  approximation: 

R/d 


P 


0.774 


(x/ d) 


2/3 


2/3 


0.121  M= 


[c„/(x/d)]“'’ 


Second  approximation  (simplified  expressions): 


R/d 

P/Pc« 


0.774 


Cp=/= /(x/d)=/'  -  1.09 

0.121  +  0.56 

(x/d)V’ 


(29) 

(33) 


(61) 

(58) 
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Axisymmetric  Flow 
First  approximation; 


p/p=0 


0.067 


Second  approximation: 
Exact; 

=  0.795  / 


+  3.15 


x/d  \ 

<  4^/ 


Simplified; 

=  0.795/[m:,  ^/^/U/d) 

“oc 

Simplified; 

p/p^  =  0.067  +  0.44 


(37) 

(41) 


(63) 


(73) 

(70) 


OTHER  BLAST  AND  STEADY-FLOW  SOLUTIONS 


The  original  Taylor's  (Ref.  3)  solution  (first  approximation)  of  a 
spherical  blast  was  extended  by  Lin  (Ref.  6)  to  the  case  of  a  cylindrical 
blast.  Lin  also  obtained  the  axisymmetric,  steady-flow  solution  (first 
approximation)  analogous  to  the  cylindrical  shock.  At  about  the  same 
time,  Sakurai  (Refs.  4,  5)  published  first  and  second  approximation 
solutions  of  plane,  cylindrical,  and  spherical  blasts.  Cheng  and  Pallone 
(Ref.  7)  applied  the  plane  wave,  first  approximation  solution  to  the  case 
of  steady  flow  past  a  blunt  plate,  for  y  =  1.4  and  1.667  .  They  also  sug¬ 
gested  that,  in  order  to  extend  the  validity  of  the  solution  to  higher 
x/d  values,  the  pressure  should  be  interpreted  as  an  increment  over  the 
free-stream  pressure  rather  than  as  the  absolute  value. 

In  Table  8,  the  numerical  values  of  coefficients  quoted  by  the  above- 
mentioned  authors  are  compared  with  the  values  given  by  Sakurai  (Ref.  4) 
and  used  in  this  paper. 
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TABLE  8 

VALUES  OF  COEFFICIENTS 

IN  THE  FIRST  BLAST-WAVE  APPROXIMATION  (y  =  1.4) 


Coefficient 

J 

0 

(o) 

g(o) 

a 

0 

1 

0 

1 

Sakurai  (Ref.  4) 

1.  696 

0.  877 

0.  455 

0.424 

Cheng  and  Pallone  (Ref.  7) 

1. 121* 

0.  319 

Lin  (Ref.  6) 

0.  858 

0.4317 

*Cheng  et  al  (Ref.  8)  pointed  out  that  this  value  is  in  error. 


Lees  and  Kubota  (Refs.  9  and  10)  applied  Sakurai's  results  to  obtain 
second  approximation,  steady-flow  solutions  analogous  to  the  case  of 
cylindrical  blast.  They  derived  "simplified"  expressions  for  the  case 
of  a  hemisphere -cylinder,  as  follows: 

p/Poo  =  0-0655  M^/(x/d)  +  0.405  (74) 

R/d  =  0.78  +  1.62  (x/d)  /  M^]  (75) 

These  expressions  correspond  to  the  binomial  expansion  of  Eqs.  (66) 
and  (63),  using  the  modified -Newtonian  drag  coefficient  for  a  hemi¬ 
sphere  at  =  7.7,  y  =  1.4  (Cj3  =  0.914). 

Love  (Ref.  11)  proposed  a  method  of  calculation  of  bluntness -induced, 
inviscid,  steady-flow  pressures  downstream  of  hemispherical  and  hemi- 
cylindrical  noses.  In  his  method,  a  value  of  the  pressure  at  the  shoulder 
is  calculated  or  assumed,  and  the  pressure  decay  as  given  by  the  first 
approximation,  blast-wave  analogy  is  taken.  Love's  expression  in  terms 
of  free-stream  pressure,  ,  is  as  follows: 

p/p.  -  [l  .  (x'/dlT'  +  [l  +  (x7d)“']^‘  <76) 

where 

b  =  Psh/Poo  =  (Psh/Po')  (Po'/Poc) 

and 

=  shoulder  pressure,  i.  e. ,  pressure  at  the  juncture  of  hemi¬ 
spherical  (or  hemicylindrical)  nose  to  the  cylindrical  (or 
slab)  afterbody; 

P„'  =  stagnation  point  pressure; 

x7d  =  dimensionless  distance  measured  from  the  shoulder. 
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Theoretical  calculations  indicate  that,  at  Mach  numbers  >  5,  Psh/Po' 
mains  approximately  constant  and  equal  to  0.  045  in  axisymmetric  flow 
and  to  0.  117  in  two-dimensional  flow  (y  =  1.4).  The  values  of  exponent  a 
are  taken  as  1  and  2/3,  respectively,  in  these  two  cases,  following  the 
blast  analogy.  Using  these  values.  Love’s  expressions  can  be  written 
as  follows : 

1.  hemicylinder  -  flat  plate: 

p/p^  =  0.117  [l  +  (x7d)^/^]“‘  (p//p^)+[l  +  (x'/d)-^/^]  "o 

(77) 


2.  hemisphere -cylinder: 

r  ,  ...  1  -  0.045  p„  7p, 

p/Poo  =  (1  +  x7d)“"Pj,Vp^  +  [l  +  (x'/d)  J 


(78) 


The  value  of  p^  7p^  is  a  function  of  the  free-stream  Mach  number. 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
WITH  APPROXIMATE  SOLUTIONS 


The  various  blast  analogy  and  other  approximations  listed  above  are 
compared  here  with  three  kinds  of  data; 

1.  results  of  theoretical  calculations  of  flow  about  simple  bodies, 
for  ideal  and  real  air; 

2.  measurements  from  wind  and  shock  tunnel  experiments; 

3.  results  of  correlations,  in  terms  of  blast-wave  analogy  param¬ 
eters,  of  theoretically  calculated  flows. 

As  regards  experimental  data,  only  those  measurements  made  at 
Reynolds  numbers  high  enough  to  render  the  viscous  effects  negligible 
are  included. 


PLANE  FLOW 

On  the  basis  of  calculation  (method  of  characteristics)  and  correla¬ 
tion  of  flow  about  flat  plates  with  sonic  wedge  leading  edges,  Baradell 
and  Bertram  (Ref.  12)  suggested  the  following  expression*  for  pressure 
distribution  (y  =  1.4): 

p/p^  =  0.117  [VI^[cD/(x/d)]^/’  +  0.732  (79) 


*This  is  derived  from  Eq.  (7)  of  Ref.  12,  with  (p  -  p^)/p^  increased 
by  3  percent  for  y  =  1.4,  as  suggested  by  the  authors. 
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This  formula  is  quite  similar  (see  Figs.  2  and  8)  to  the  simplified  sec¬ 
ond  approximation,  Eq.  (58),  and  gives  slightly  larger  values  of  p/p^  at 

small  values  of  [Cp/ (x/d)]^-^’,  and  vice  versa.  The  formula  was  de¬ 
rived  on  the  basis  of  two  calculations,  both  for  «  1.4  and  for  =  9.5 
and  =  20. 

Pressure  distribution  measured  on  a  flat  plate  with  a  semi -cylindrical 
leading  edge  is  compared  with  various  approximations  in  Fig.  2.  The 
experimental  data  extend  to  x/d  =  10.5.  Within  this  range,  all  blast-wave 
approximations  calculated  for  the  experimental  pressure  drag  coefficient  of 
Cp  =  1.27  show  reasonable  agreement  with  experimental  values  at  x/d  >  1.5. 
The  correlation  formula  of  Ref.  12,  Eq.  (79),  gives  values  very  slightly 
larger  than  the  second  approximation.  Although  Love's  method  predicts 
correctly  the  shoulder  pressure,  it  is  least  accurate  at  intermediate 
x/d  values . 

Also  in  Fig.  2  a  curve  is  included  showing  second  blast  approxima¬ 
tion  computed  for  =  8  and  =  1  instead  of  the  experimental  value 
of  1.  27.  The  agreement  with  measurements,  although  presumably  for¬ 
tuitous,  is  excellent. 

In  Figs.  3,  4,  and  5,  approximations  to  the  pressure  distribution 
and  shock  shape  in  plane  flow  are  compared  with  the  theoretical  results 
obtained  from  calculations  using  the  method  of  characteristics. 

The  superiority  of  the  second  approximation  to  the  pressure  distri¬ 
bution  at  large  x/d  values  for  a  configuration  having  a  high  nose  drag  is 
very  evident  at  =  6.86  from  Figs.  3a  and  b.  The  difference  between 
the  theoretical  and  approximate  pressure  ratios  at  any  given  x/d  is  sub¬ 
stantially  constant,  and  the  approximate  pressure  ratio  is  within  20  per¬ 
cent  of  the  theoretical  value  even  at  x/d  =  150.  A  still  better  agreement, 
to  within  10  percent,  is  obtained  with  the  correlation  formula,  Eq.  (79). 

In  Fig.  4,  the  corresponding  theoretical  and  approximate  shock 
shapes  are  shown,  and  the  superiority  of  the  second  blast  approximation 
is  again  evident.  This  is  also  the  case  at  =  12.3,  Fig.  5;  the  second 
blast  approximation  gives  values  to  better  than  10  percent  of  the  experi¬ 
mental  observations. 

Results  of  flow  field  computations  for  real  air  obtained  by  the  General 
Electric  Company  (Refs.  16  and  17)  are  shown  in  Figs.  6  and  7.  The 
blast  analogy  does  not  predict  the  pressure  distribution  accurately  for 
y  =  1.4;  a  much  better  agreement  was  obtained  for  y  =  1.3,  Fig.  6.  The 
shock  slope,  as  in  the  previous  cases,  is  well  approximated  by  blast 
analogy  for  y  =  1.4.  In  Fig.  7,  the  shock  Mach  numbers,  Mjj,  computed 
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from  Eqs.  (24)  and  (31),  are  shown  and  indicate  that  real  gas  effects 
might  affect  a  significant  portion  of  the  flow  field. 

The  data  shown  in  Figs.  2  through  7  are  plotted  in  terms  of  blast 
analogy  correlation  parameters  in  Figs.  8  and  9. 

All  theoretical  results  for  pressure  distribution,  for  y  =  1.4,  Fig.  8, 
correlate  well  with  Eq.  (79).  The  experimental  data  fall  below  the  cor¬ 
relation  of  theoretical  data.  The  inadequacy  of  the  first  blast  approxi¬ 
mation  is  again  clearly  seen. 

Figure  8  indicates  that,  although  blast  analogy  provides  suitable 
correlation  parameters,  the  range  of  its  applicability,  in  terms  of  x/d, 
depends  on  values  of  Cp  and  M^.  For  example,  with  Cp  =  1.4,  the  pre¬ 
diction  of  pressure  ratio  becomes  poor  at  =  9.5  when  (x/d)^/V(M^  Cp 
<  0.025  (x/d  <  5),  whereas  at  =  20,  it  is  good  down  to  about  0.004(x/d  =  3). 

Comparison  of  results  of  theoretical  calculations  for  ideal  (y  =  1.4) 
and  real  air,  for  Mach  numbers  of  20  and  15,  respectively,  presumably 
indicates  the  already  noted  significant  real  gas  effects. 

In  Fig.  9,  the  available  data  on  shock  shape  are  plotted  in  terms  of 
correlation  parameters  of  the  second  approximation  blast  analogy,  to¬ 
gether  with  the  theoretical  curve.  The  two  sets  of  data  correlate  poorly. 
The  experimental  results  agree  well  with  the  theoretical  curve  at  large 
values  of  x/d. 

From  the  above  review  of  the  limited  data  available  for  plane  flow, 
it  is  apparent  that  blast  analogy  provides  useful  correlation  of  pressure 
distributions  and  shock  shapes.  Pressures  are  well  predicted  by  the 
second  blast  analogy  approximation  for  ideal  gas,  y  =  1.4;  results  of  a 
calculation  for  real  air  at  =  15  indicate  pressures  lower  by  20  per¬ 
cent  or  more.  Shock  slopes  are  well  predicted  by  the  second  blast  analogy 
approximation,  whereas  the  shock  location  shows  an  outward  displacement 
relative  to  blast  analogy  prediction. 


AXISYMMETRIC  FLOW 


Experimental  Data 


The  available  experimental  pressure  distributions  on  hemisphere 
cylinders  at  Mach  numbers  6,  7.  7,  and  8  are  compared  with  blast  analogy 
and  Love's  (Ref.  11)  predictions  in  Figs.  10a,  b,  and  c.  In  all  cases  at 
x/d  >  3,  the  second  blast  analogy  approximation  gives  values  much  more 
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realistic  than  the  first  approximation  and  predicts  the  pressure  distribu¬ 
tion  accurately  at  M^=  7.7  and  8  .  Love's  method  gives  the  pressure  level 
accurately  at*M^  =  6,  x/d  >  4,  but  is  less  accurate  than  the  second  blast 
analogy  approximation  at  =  8. 

Pressure  distribution  on  a  flat-nose  cylinder  at  =  8  is  shown  in 
Fig.  11.  Within  the  limited  range  of  measurements  (x/d  <  5.5),  both  first 
and  second  blast  analogy  approximations  give  an  accurate  estimate  of  the 
pressure  on  this  high  drag  (Cp  =  1.64)  body. 

The  experimental  and  blast  analogy  data  on  shock  shapes  are  given 
in  Fig.  12  for  two  of  the  cases  already  considered.  As  in  the  plane  flow, 
the  second  approximation  predicts  closely  the  shock  slope. 

Lees  and  Kubota  (Refs.  9  and  10)  measured  the  radial  pitot  pres¬ 
sure  distribution  around  the  hemisphere-cylinder  model  at  =  7.7  at 
station  x/d  =  3  and  computed  the  corresponding  distributions  of  other 
quahtities.  These  results  for  p/p^  and  p/p^  are  shown  in  Figs.  13  and  14 
together  with  the  corresponding  blast  analogy  predictions.  The  latter  were 
obtained  from  Eqs.  (3)  and  (8)  with  y  =  l/M^  computed  from  Eqs.  (24),  (39), 
and  (64). 

The  experimental  pressure  distribution.  Fig.  13,  shows  an  excellent 
agreement  with  the  second  blast  approximation  over  the  whole  range  of 
r/R.  At  the  shock,  the  observed  pressure  ratio  is  5,  whereas  the  one 
computed  from  the  second  approximation  blast  analogy  is  5.13. 

The  density  distributions  are  compared  in  Fig.  14.  As  mentioned 
previously,  the  second  blast  approximation  gives  the  exact  density  ratio 
across  the  shock  only  at  high  shock  Mach  numbers.  In  this  case,  Mj^  =  2.13, 
and  hence  a  rather  poor  agreement  is  obtained  with  the  experimental  dis¬ 
tribution.  However,  the  exact  p/p^  ratio  for  Mj^  =  2.13  equals  2.  85  and 
agrees  well  with  the  experimental  value  (Fig.  14). 

Nevertheless,  in  view  of  the  good  agreement  in  the  pressure  distri¬ 
butions  in  Fig.  13,  it  can  be  considered  that  the  flow  field  closely  approxi¬ 
mates  the  second  approximation  blast  solution.  It  is  then  possible  to  write, 
for  the  density  distribution,  an  expression  which  is  exact  at  the  shock  and 
to  determine  a  function  x'(«:)  from  the  experimental  data.  Thus,  by  analogy 
to  Eq.  (17), 

p/Poo  =  h(«:)  /[i  +  x'Wy]  (80) 

(o) 

where  h(/c)  is  given  by  the  first  blast  approximation  and  x'  is  computed  to 
fit  the  experimental  density  distribution  of  Fig.  14.  The  values  of  x'  are 
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given  in  Table  9  and  were  computed  assuming  Mj^  =  2.11,  which  corre¬ 
sponds  to  the  experimental  pressure  and  density  ratios  across  the  shock. 
At  higher  shock  Mach  numbers,  Eqs.  (80)  and  (8)  are  in  close  agreement. 


TABLE  9 

VALUES  OF  x'  (Mj^  =  2.11) 


r/R 

1 

0.  95 

o 

CD 

O 

00 

0.  7 

CD 

O 

0.  5 

0.4 

x' 

5 

1.  75 

0.  1 

-1.  45 

-2.  45 

-3.  2 

-3.  8 

-4.  25 

Theoretical  Data  for  y  =  1.4 

Van  Rise  (Ref.  14)  published  results  of  calculations  of  flow  fields 
about  cylinders  with  a  variety  of  nose  shapes,  ranging  in  drag  coefficients 
from  0.  037  to  1.  37.  He  used  the  method  of  characteristics  to  obtain  sur¬ 
face  pressure  distributions  and  shock  shapes  in  the  flow  of  perfect  gas 
with  y  =  1.4  and  5/3.  In  most  cases,  the  nose  drag  coefficient  was  com¬ 
puted  from  the  surface  pressure  distribution*. 

Pressure  distributions  obtained  by  Van  Rise  are  shown  in  Figs.  15a,  b, 
and  c  in  terms  of  the  blast  analogy  correlation  parameter,  (x/d)/(M:,  ^). 
Also  in  these  figures  is  included  the  curve  given  by 


p/p^  =  0.06  /(x/d)  +  0.55  (81) 

which  was  derived  by  Van  Rise  as  representing  the  best  correlation  of 
all  of  his  results  for  air  (y  =  1.4).  Equation  (81)  above  closely  approxi¬ 
mates  the  theoretical  blast  analogy  expression,  Eq.  (70)  (see  Fig.  24). 

It  is  immediately  apparent  from  Van  Rise's  results  that  the  blast 
analogy  provides  an  excellent  means  of  correlating  the  pressure  data 
over  a  wide  range  of  drag  and  Mach  numbers.  It  is  also  apparent  that 
the  range  of  validity  of  an  expression  such  as  Eq.  (81)  is  limited,  in  each 
case,  to  a  certain  range  of  the  correlation  parameter.  Based  on  Figs.  15b 
and  c,  for  bodies  with  »  l  and  at  >  10,  Eq.  (81)  predicts  accurately 

the  pressure  distribution  at  x/d  >  2.  The  lower  limits  (in  terms  of  free- 
stream  Mach  number  and  correlation  parameter,  (x/d)  /  ,  of  validity 

of  blast  analogy  correlation  are  indicated  in  more  detail  in  Fig.  16  for 
different  drag  coefficients,  based  on  agreement  in  the  pressure  ratio  to 


^Values  of  Cj^  assumed  rather  than  computed  are  given  in  parentheses 
in  Figs.  15  and  18.  Only  the  results  for  y  =  1.4  are  here  considered. 
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within  10  percent.  The  theoretical  data  of  Van  Hise  is  self-consistent, 
whereas  the  experimental  data  indicate  higher  limits  of  the  correlation 
parameter  for  given  Cp  and  .  Similar  data  are  presented  in  Fig.  17 
for  =  40  and  a  wide  range  of  Cp. 

As  regards  the  upper  limits  of  the  validity  of  the  correlation  param¬ 
eter,  the  available  data  are  not  extensive  enough  to  allow  their  determina¬ 
tion.  In  the  majority  of  cases  considered,  good  correlation  was  obtained 
for 

(x/d)/(M^^)<  0.1  or  p/p^  >  1. 

In  a  number  of  cases  in  which  calculations  have  been  made  into  the  region 
of  p/p^  <  1  (Figs.  15b  and  c,  at  =  10  and  6.9),  significant  deviations 
from  Eq.  (81)  are  evident  at  p/p^  <  1,  but  the  data  appear  nevertheless 
to  be  well  correlated  by  the  blast  analogy  parameter.  This,  however,  is 
not  generally  true,  as  indicated  in  Fig.  15c  by  the  pressure  distribution 
at  =  20,  Cq  =  1.37,  which  deviates  from  other  data  already  at  p/p^  =  1.3 
and  at  smaller  p/p^  values.  As  pointed  out  by  Van  Hise*,  the  results  at 
p/p^  <  1.3  may  be  subject  to  significant  cumulative  errors. 

The  corresponding  available  results  for  the  shock  shape  are  shown 
in  Fig.  18  in  terms  of  the  blast  analogy  correlation  parameters.  The 
correlation  is  good  considering  the  wide  range  of  drag  and  Mach  numbers 
and  falls  slightly  above  the  predicted  (R/<1)/(M^  )  values. 


Theoretical  Data  for  Real  Air  in  Equilibrium 


Two  sets  of  theoretical  calculations  of  flow  fields  around  a  hemisphere- 
cylinder  in  real  air  in  equilibrium  are  available  for  comparison  with  blast 
analogy  predictions  and  correlations.  Feldman  (Ref.  15)  published  results 
of  characteristic  calculations  for  «  18  (17,  500  ft/sec)  at  60,  000-ft  alti¬ 
tude,  and  the  General  Electric  Company  (GE)  (Ref.  16)  computed  similar 
data  for  =  15  to  19  at  100,  000  to  200, 000-ft  altitude,  using  the  flow 
field  method  of  Gravalos  (Ref.  17).  Details  of  the  conditions  for  the 
GE  calculations  are  given  in  Table  10  on  the  following  page. 


^Private  communication. 
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TABLE  10 

GE  CALCULATIONS  OF  FLOW  FIELDS  AROUND  A  HEMSPHERE-CYLINDER 
IN  REAL  AIR  IN  EQUILIBRIUM 


M 

oo 

ft/sec 

h 

CX)  ’ 

ft 

Poo. 

Ib/ft2 

Poo  X  10% 

slug/ ft  ^ 

T 

(15)  ^ 

14,  532 

0.  897 

100,  000 

23. 085 

32. 114 

418. 79 

18.  24 

20,  000 

0.  914 

150,  000 

3.  0597 

3.  5642 

500.  11 

18.  1 

20,  000 

0.  91 

175,  000 

1. 2334 

1.  4123 

508. 79 

19.  25 

20,  000 

0.  908 

200,  000 

0. 47151 

0.  6118 

449.  0 

Shock  shapes  and  pressure  distributions  obtained  by  Feldman  (Ref.  15) 
and  from  blast  analogy  predictions  (for  y  =  1.4,  =  i)  are  shown  in  Fig.  19. 

As  for  the  other  data,  the  shock  shape  agrees  well  with  the  second  blast 
approximation,  except  for  a  substantially  constant,  outward  displacement 
of  about  (1/4)  d. 

The  pressure  distribution,  as  given  by  the  second  blast  analogy  ap¬ 
proximation,  shows  a  fair  agreement  with  theory  at  large  x/d  values. 

The  GE  data  are  shown  in  Figs.  20  through  23  for  =  15  and  18.1. 

The  shock  shapes  and,  more  particularly,  the  shock  slopes  are  well  pre¬ 
dicted.  At  =  15  (Fig.  21),  the  pressure  distribution  agrees  well  with 
all  blast  analogy  approximations  for  y  =  1.4  at  x/d  >3.  At  =  18.1  (Fig.  23), 
the  agreement  is  somewhat  poorer  at  x/d  <  4,  probably  indicating  a 
stronger  influence  of  real  gas  effects.  However,  compared  to  the  results 
obtained  by  Feldman  (Fig.  19b),  a  much  better  agreement  with  the  blast 
analogy  is  evident  at  small  x/d  values. 

For  comparison,  results  for  y  =  1.3  are  also  included  in  Figs.  21 
and  23.  It  is  apparent  that  a  change  in  y  does  not  account  for  the  dif¬ 
ferences  between  the  calculated  and  predicted  pressure  distributions. 

The  real  air  data  are  shown  in  terms  of  the  blast  analogy  parameters 
in  Figs.  24,  25,  and  26,  together  with  the  experimental  results  and  curves 
representing  blast  analogy  solutions  and  other  correlations. 

At  lower  values  of  the  pressure  ratio  (Fig.  24)  the  real  air  results 
show  a  fair  agreement  with  the  other  data.  At  pressure  ratios  above 
about  4,  the  real  air  results  for  =  18  to  20  show  a  marked  divergence 
from  the  blast  analogy  and  the  correlation  of  theoretical  data  for  y  =  1.4. 
Also,  unlike  the  theoretical  computations  for  cone-cylinders  and  cone- 
sphere -cylinders  (Fig.  15),  the  pressure  distributions  at  =  18  to  20 
in  real  air  correlate  well  from  the  shoulder  (x/d  =  0.5)  downstream. 
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These  trends  are  more  clearly  shown  in  Fig.  25.  It  is  evident  that, 
whereas  the  =  15  results  in  real  air  are  in  fair  agreement  with  the 
y  =  1.4  characteristic  and  blast  analogy  correlations,  the  =  18  to  20 
data  correlation  shows  a  large  change  in  slope,  which  cannot  result 
from  a  simple  change  in  the  value  of  y. 

Shock  shapes  for  real  air  are  compared  in  Fig.  26  with  the  theoreti¬ 
cal  blast  analogy  approximations  (y  =  1.4)  and  experimental  data  (y  =  1.4). 

At  larger  values  of  x/d  [(x/d)  /(m^  )  >  O.Ol],  the  real  air  results  closely 

approach  the  second  approximation  blast  analogy  prediction  for  y  =  1,4. 

As  indicated  in  Fig.  26,  the  agreement  would  be  worse  for  y  <  1.4.  The 
only  available  experimental  data  fall  above  the  blast  analogy  prediction, 
as  already  indicated  (Fig.  12b).  At  small  x/d  values,  smaller  R/d  values 
and  smaller  shock  slope  than  predicted  are  found.  This  trend  for  real 
air  coincides  with  results  of  the  characteristic  computations  for  y  =  1,4 
(Fig.  18). 

The  radial  pressure  and  density  distributions  computed  by  GE  (Ref.  16) 
at  =  18.1  are  compared  with  blast  analogy  predictions  for  y  =  1.4  in 
Figs.  27  and  28. 

As  already  observed,  the  pressure  distribution  (Fig.  27)  is  well  pre¬ 
dicted  at  larger  x/d  values;  in  all  cases,  the  first  approximation  agrees 
better  with  the  theoretical  results  than  the  second  approximation.  This 
would  be  expected  in  view  of  the  similarly  better  agreement  in  shock 
slope  (Fig.  22). 

The  density  distribution  (Fig.  28)  is  in  all  cases  well  predicted  at 
r/R  >  0.8.  At  the  shock,  the  second  approximation  gives  more  accurate 
values,  and  still  closer  ones  are  obtained  using  Eq.  (80).  At  r/R  <  0.8, 
the  predicted  density  is  much  smaller  than  the  computed  one. 


CONCLUDING  REMARKS 


Comparison  of  blast  analogy  (second  approximation)  solutions  with 
correlations  based  on  characteristic  calculations  for  plane  and  axisym- 
metric  flows  shows  good  agreement  between  the  two  and  thus  provides  a 
confirmation  of  the  validity  of  the  blast  analogy  method.  The  range  of 
applicability  of  the  blast  analogy  solutions  depends  on  the  values  of  param¬ 
eters  (e.  g. ,  and  Cq  )  which  define  the  flow  field;  however,  in  general, 
at  high  Mach  numbers  and  for  ideal  gas  (y  =  1.4),  the  blast  analogy  pre¬ 
dicts  rather  accurately  the  pressure  distributions  and  shock  shapes  from 
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2  to  3  diameters  downstream  of  the  nose  to  downstream  stations  at  which 
the  pressure  has  decayed  to  the  free-stream  value.  Although  the  shock 
is  displaced  from  the  body  more  than  predicted  by  blast  analogy,  the 
latter  gives  accurate  values  of  shock  slope. 

As  regards  application  of  the  blast  analogy  method  to  flow  of  real  air, 
the  available  data  are  limited  to  semi-cylinder-plate  and  hemisphere- 
cylinder  configurations  at  =  15  to  19  ,  at  altitudes  from  60,  000  to 
200,  000  ft.  Comparisons  indicate  that  blast  analogy  for  y  =  1.4  predicts 
accurately  the  shock  slope;  the  pressure  distribution  is  predicted  about 
as  well  as  for  experimental  measurements  made  at  lower  Mach  num¬ 
bers  (6  to  8)  and  temperatures. 
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•  VKF-AEDC,  flat  plate  with  a 
semi-cylindrlcal  leading  edge 
Re^  -  0.17  X  10^,  air,  =  8.08 

Cj)  °  1,27,  fron  nose  pressure  dlstributioD 
Cp  °  1.22,  modified  Newtonian 

—  —  —  —  let  approx,  Eq.  (33)  j  y  -  1.4 

— 2nd  approx,  Eq,  (58)  |  =  1.27 

+  Eq.  (79)  (Ref.  12)  I  -  8.08 

— Love's  method,  Eq,  (77),  -  8.08,  Ref.  11 

—  ^  2nd  approx,  Eq.  (58).  7  =  1.4,  (V,  ■=  1. 


2nd  approx,  Eq.  (58),  y 
M  «  8 
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Bertram  and  Baradell,  Fig.  2,  Ref.  13, 
flat  plate  with  a  sonic  wedge  (43~deg 
haif- angle)  leading  edge;  characteristic 
solution,  3rd  approx, 
r  “1.4,  Cp-  1.374,  M^-6.86 

2nd  approx,  V  v  =  1. 4 
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a.  x/d  <  15 

Fig.  5  Experinwntal  and  Predicted  Shock  Shapes 
in  Plane  Flow,  =  12.3 
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Cheng,  et  al,  Fig.  4,4, 
Ref.  8,  flat  plate  with 
a  square  leading  edge 
711  <  Re^i  <  15,000 
0.0096  in.  <  d  <  0.203  in. 
M  =12.3,  air 


10  20  30 

x/d 


b.  x/d  <  dO 

Fig.  5  Continued 


Cheng,  et  al,  Fig.  4.4, 
Ref.  8,  flat  plate  with 
a  square  leading  edge 
711  <  Red  <  15,000 
0.0096  in.  <  d  <  0.203  in 
M  «  12.3,  air 

2nd  approx,\  y  =  1.4 
Cb  =  2 


1st  approx, 
Eq.  (29) 


M  =  12.3 

00 


a 

n 


2  3  4  5  6  7 

x/d 


Fig.  7  Shock  Shopes  In  Plono  Flow  Calculated  by  ihe  Flow  Field  Method  for  Real  Air  and  by  Blast  Analogy,  =  15 
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FiQ*  9  Shock  Shopo  in  Plano  Flow  Corrolatod  by  Blast  Analogy  Paramotors 
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#  Ues  and  Kubota,  n 

Refs.  9  and  10,  li 

hemisphere  cylinder,  ^ 

Red  ~  d- 13  X  10^,  air,  ^ 

Mco“7.7.  Z 


Cp"  0.914^  mod.  Newt,  y-1.4 
—  —  2nd  approx,  \  y-l.4 

CD-a914 

- 1st  approx,  „  -- 

Eq.  Ml) 


Fig.  10 


•  VKF-AEDC,  hemisphere  cylinder,  air. 

Red  ->  1.34  X  10^,  Ma>  -  8l08 

Cq  *  Oi09,  from  pressure  distribution 

Cd  -  a  914,  modified  Newtonian  (r  =  1.4) 


2nd  approx,  V  y  >=■  1.4 
CD-a» 

S!«r  ) 


J - 1 - 1 _ 

12  3 

x/d 

11  Expvrimenial  aii<l  Pradielwd  Pressure  Distribution! 


I  VKF-AEDC,  flat-nose  cylinder 

(Symbols  indicate  spread  of  pres-| 
sure  measured  at  45-deg  intervals 
of  circumference.) 

Re^  =  1,01  X  10«>,  Air,  H  =8.08 

00 

Cj)  =  1.64,  from  pressure 
distribution 

Cd  =  1.83,  modif.  Newtonian 


2nd  approx, 
Eq.  (70) 


_  1st 
Eq. 


approx 

(41) 


A  y=  1.4 

Ud  =  1.6 

’)  M  =  8.0 
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08 


I 
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Cylinder,  =  8 
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p/p« 


15  Correlation  of  Theoretical  Pressure  Distribution  in  Ax i symmetric  Flow  for  y  =  1-4,  Ref.  14 
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b.  Cd  0.79  an 
Fig.  15  Contin 


n  m  R7I  \  Cone-Hemisphere-Cylinder 
A  (0  87)  10  )  530  Half-Vertex  Angle 


•  0  79  40  ^ 

■  rn  70i  iR  (  Cone-Hemisphere-Cylinder 

A  (o!  79)  10  )  Half-Vertex  Angle 
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0.0001 
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Fig.  16  Limift  of  Validity 
y  **  1 .4  (  Based  or 


Fig.  17  Limit  of  Validity  of  Blost  Analogy  Correlation  at  -  40,  y  =  1.4 
(Based  on  Computations  by  Van  Hise,  Ref.  14) 


Fig.  18  Correlation  of  Theoretical  (Method  of  Characteristics,  y  «  1.4)  Shock  Shapes  in  Axisymmetric  Flow 
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Feldman,  Ref.  15, 
hemisphere-cylinder , 
method  of  character¬ 
istics,  real  air, 
60,000  ft  altitude, 
17,500  ft/sec,  ~ 

2nd  approx,  \  7  •=  1- 

Eq.  (63)  (cp  -  1 

1st  approx,  I  „  , , 

Eq.  (37)  /“oo 


a.  Shock  Shapes 

Fig.  19  Shock  Shopos  and  Prossuro  Disiributiens  Calculated  by  the  Method  of  Characteristics  and 
Blast  Analogy  in  Axi symmetric  Flow,  ••  18 


x/d 

b.  Pressure  Distributions, 


Fig.  19  Continue 


Feldman,  1960, 
hemisphere-cylinder , 
method  of  character¬ 
istics,  real  air, 
60,000  ft  altitude, 

17,500  ft/sec,  M  “  18 

00 

2nd  approx, 

Eq.  (70) 

1st  approx, 

Eq.  (41) 


00 


■■■■■■  Real  Air,  Ref.  16 

1st  approx,  ' 

Eq.  (41) 

■  2ad  approx, 

Eq.  (70)  “  “ 

van  HisCj  I 

Eq.  (81) 

■  lat  approx,  y  - 
Kq.  (40) 
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Fig,  21  Pressure  Distribution  on  a  Hemisphere  Cylinder  at  ^15  in  Real  Air 


1st  approx, 
Eq.  (37) 


Fig«  23  Preisure  Diitribution  on  a  Hemisphere  Cylinder  at  18.1  in  Real  Air 
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Fig.  24  Correlation  of  Experimental  and  Theoretical  (for  Real  Air)  Pressure  Distribution  in  Axisymmetric  Flow 
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FTg,  25  Corralation  of  Pressure  Disfribution  In  Axisymmetric  Flow  of  Reol  and  Ideol  [y  -■  1.4)  Air  ^ 
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Radial  Prassura  Disiributioii  in  Flow  of  Real  Air 


